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Abstract We present a nonlinear model of weakly curved rod, namely the 
type of curved rod where the curvature is of the order of the diameter of the 
cross-section. We use the approach analogous to the one for rods and curved 
rods and start from the strain energy functional of three dimensional nonlinear 
elasticity and do not presuppose any constitutional behavior. To derive the 
model, by means of .T-convergence, we need to propose how is the order of 
strain energy related to the thickness of the body h. We analyze the situation 
when the strain energy (divided by the order of volume) is of the order /i^. 
That is the same approach as the one when Foppl-von Karman model for 
plates and the analogous model for rods are obtained. The obtained model is 
analogous to Marguerre-von Karman for shallow shells and its linearization is 
the linear shallow arch model which can be found in the literature. 

Keywords weaky curved rod ■ Gamma convergence ■ shallow arch • 
asymptotic analysis 
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1 Introduction 

The study of thin structures is the subject of numerous works in the theory 
of elasticity. There is a vast literature on the subject of rods, plates and shells 
(see [51I81I9]). 

The derivation and justification of the lower dimensional models, equilib- 
rium and dynamic, of rods, curved rods, weakly curved rods, plates and shells 
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in linearized elasticity, by using formal asymptotic expansion, is well estab- 
lished (see [SJH] and the references therein) . In all these approaches one starts 
from the equations of three-dimensional linearized elasticity and then via for- 
mal asymptotic expansion justify the lower dimensional models. One can also 
obtain the convergence results. In [2111] the linear model of weakly curved rod 
(or as it is called shallow arch) is derived and the convergence result is ob- 
tained. We call weakly curved rods or shallow arches those characterized by 
the fact that the curvature of their centerline should has the same order of 
magnitude as the diameter of the cross section, both being much smaller than 
their length. 

Formal asymptotic expansion is also applied to derive non linear models of 
rods, plates and shells (see j8ki9u22j and the references therein), starting from 
three-dimensional isotropic elasticity (usually Saint- Venant-Kirchoff material). 
Hierarchy of the models is obtained, depending on the the order of the external 
loads related to the thickness of the body h (see also [Tl] for plates). 

However, formal asymptotic expansion does not provide us a convergence 
result. The first convergence result, in deriving lower dimensional models from 
three-dimensional non linear elasticity, is obtained applying /^-convergence, 
very powerful tool introduced by Degiorgi (see [51ITU]). Using r'-convergence, 
elastic string models, membrane plate and membrane shell models are obtained 
(see pl lTBlfTT] ). It is assumed that the external loads are of order The ob- 
tained models are different from those ones obtained by the formal asymptotic 
expansion in the sense that additional relaxation of the energy functional is 
done. 

Recently, hierarchy of models of rods, curved rods, plates and shells is 
obtained via F- convergence (see [ni[n|Pll|TCl[Tgi[^[^|5ni|5T] ) . Influence of 
the boundary conditions and the order and the type of the external loads is 
largely discussed for plates (see [T21[TS] 1. Let us mention that T-convergence 
results provide us the convergence of the global minimizers of the total energy 
functional. Recently, compensated compactness arguments are used to obtain 
the convergence of the stationary points of the energy functional (see [^[25] ). 

Here we apply the tools developed for rods, plates and shells to obtain 
weakly curved rod model by /^-convergence. It is assumed that we have free 
boundary conditions and that the strain energy (divided by the order of vol- 
ume) is of the order /i"*, where h is the thickness of the rod. This corresponds 
to the situation when external transversal dead loads are of order (see Ke- 
mark|5]). The order /i* of the strain energy gives Foppl-von Karman model for 
plates, Marguerre-von Karman model for shallow shells the analogous model 
for rods (see |14[[^I5^ ). The obtained model is non linear model of the lowest 
order in the hierarchy of models and its linearization is shallow arch model, 
obtained in for isotropic, homogenous case (see for comparison Remark 
[7]d)). Here we do not presuppose any constitutional behavior and thus work 
in a more general framework. The main result is stated in Theorem [5] 

Throughout the paper A or {A}~ denotes the closure of the set. By a 
domain we call a bounded open set with Lipschitz boundary. I denotes the 
identity matrix, by SO (3) we denote the rotations in M.^, by so(3) the set of 
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antisymmetric matrices 3x3 and M.^y^ denotes the set of symmetric matrices. 
By sym A we denote the symmetric part of the matrix, sym A = ^(A + A"^). 
ei, 62, 63 are the vectors of the canonical base in R'^. By Vh we denote V/,. = 
Vei + -^Vea.e.,- ||/||ci(f2) stands for norm of the function / : i7 C M" — > 
M i.e. \\f\\c^^2) = max^jgr? |/| + J27=i'^^^xen \dif\. -> denotes the strong 
convergence and ^ the weak convergence. 



2 Setting up the problem 

Let C be an open set having area equal to A and Lipschitz boundary. 
For all h such that < ft. < 1 and for given L we define 

uj^ = hoj, Q'^ = (0,i) X huj. (2.1) 

We shall leave out superscript when /i = 1, i.e. Q = , ljj ~ uj^ . Let us by 
/j,(a;) denote 

/i(w) = / {x\ + x\)dx2dxj,. (2-2) 



Let us choose coordinate axis such that 

/ X2dX2dX3 — / X3dX2dX3 = / X2X3dX2dX3 — 0. (2-3) 

J UJ J uj J U) 

For every h we define the curve of the form 

= {e'^ixi) = {xi,el^{xi),e'i{xi)) ew^-.x^^ (o,l)}. (2.4) 

where 9^{xi), for A; = 2,3, are given functions satisfying e C'^(0, L). Let 
{t^ , n^^b^) be the Frenet trihedron associated with the curve C'' 

i'^ = ^ (1, (e'i)', ie'})'), (2.5) 

b'' ^t'' X n'\ (2.7) 

We suppose n'* e C^{0,L) which is satisfied if {O2)" , (6*3)" do not vanish at 
the same time (which is equivalent to the fact that the curvature of is 
strictly positive for any xi £ {0,L)). The case where has null curvature 
points can be treated in the same fashion, provided that we suppose that along 
these points we have the same degree of smoothness as before with t'^, n'* and 
b appropriately chosen (see Remark [3]). We define the map 
&^{fi^') = C R^ where :== &^{[2^), in the following manner: 

©''(x'*) = (xi,02^a;i),0j(a;i))+xV(^i) + 4b''(^i) (2-8) 
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and we assume that 0^ is a diffeomorphism which can be proved if h is 
small enough and 6'^, for k = 2,3, are of the form considered here. Namely, we 
take 02 = h92, 6*3 = hO^, where Ok G C^(0,L). Let us suppose 

{{e,)"f{x,) + {{62)"f{x,)^Q, (2.9) 

for all xi G (0,-L). A generic point in or {Q^}^ will be denoted by x^ — 

{xi, X2,X^). 

Like in |241[^[T^[T51[H1[T51[T^ we start from three dimensional non linear 
elasticity functional of strain energy (see [7] for an introduction to non linear 
elasticity) 

^'(y)-=r2/' W\x'\Vy)dx\ (2.10) 

It is natural to divide the strain energy with h? , since the volume is vanishing 
with the order of h? . We are interested in finding I^-limit (in some sense i.e. 
in characterizing the limits of minimizers) of the functionals -^I^ ■ The reason 
why we divide with /i"* is that we want to obtain theory analogous to Foppl-von 
Karman for plates and rods (see |13l[T4l[25] ) and Marguerre-von Karman for 
shallow shells (see [35]). We do not look the total energy functional because the 
part with the strain energy contains the highest order derivatives (at least for 
the external dead loads) and thus makes the most difficult part of the analysis 
(see Remarks [5] and ITUl) . We shall not impose Dirichlet boundary condition 
and assume that the body is free at the boundary. The consideration of the 
other boundary conditions is also possible. We rewrite the functional on 
the domain fl, i.e. we conclude 

I^{y):= f W^'{&''oP''\x),{Vy)o0'''oP^)Aci{{V0^)oP''\x))dx, (2.11) 
J n 

where by P'' : M'^ — > M'' we denote the mapping P^{xi,X2, X3) — {xi, hx2, hxs). 
{Vy) o 0^ o P'' denotes Vy evaluated at the point 0^{P^{x)). We assume 
that for each h it is valid 

det((V6>'') oP''(a;))P^''(0''- oP''(a;),F) W{x,F),\/x S VF G R3^^ 

(2.12) 

where the stored energy function W is independent of h and satisfies the 
following assumptions (the same ones as in |25j): 

i) W : f2 X M."^^^ [0, +00] is a Caratheodory function; for some 6 > the 
function F t-^ W{x,F) is of class for dist(F, S0(3)) < 6 and for a.e. 

X e Q; 

ii) the second derivative is a Caratheodory function on the set J7 x {F G 
]^3x3 . (jist(F, S0(3)) < 6} and there exists a constant 7 > such that 

< 7|G|2 if dist(F,S0(3)) < 5 and G G M.^,^^; 



5F2 



{x,¥)[G,G] 
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iii) W is frame-indifferent, i.e. W{x,F) = W{x,'R,F) for a.e. a; G J7 and every 
FGR^^'\ReS0(3); 

iv) W{x,F) = if F e S0(3); W{x,F) > C dist^(F, S0(3)) for every F e 
M3x3^ where the constant C > is independent of x. 

Under these assumptions we first show the compactness resuh (Theorem [3]) 
i.e. we take the sequence y'* € M^^'^ M'') such that 

hmsupft, — > ^j^l'^ < +00 

and conclude how that fact affects the hmit displacement. In Lemma [2] we 
prove the lower bound, in Theorem 2] we prove the upper bound and that 
enables us to identify limit functional (Theorem [SJ . First we start with some 
basic properties of the mappings 0^ which are necessary for further analysis. 



3 Properties of the mappings 0^ 

We introduce for k — 2,3, 

Pkixi) = ^fc(^i) _ (3.1) 

Notice that 

P2+P3 = Ij = 0- (3-2) 

Let us denote p — P2P3 — P^P^ ■ 
Theorem 1 Let the functions 9^ be such that 

e'l^ixi) = hOkixi), for aU xi G (0,L),fc = 2,3. 

where 9k G C^{Q,L) is independent of h. Then there exists Hq — hQ{9) > 
such that the Jacobian matrix ^0'^{x^), where the mappings 0^ are defined 
with 112.8]) . is invertible for all x^ G fi^ and all h < ho. Also there exists C > 
such that for h < ho we have 

detV©'* = 1 + M'^(a;'*), (3.3) 



t^'ixi) ^ei + h9'^(xi)e2 + h9'^{xi)e:i + h^Oi{xi), (3.4) 
n^{xi) = P2{xi)e2 + P3{xi)ej, - h{9'2P2 + 9'^p3){xi)ei 

+h^02{xi), (3.5) 
b''{xi) = -~P3{xi)e2 +P2{xi)e3 + h{9'2P3 - 9'^p2){xi)ei 

+h^03{xi), (3.6) 
W0''{x'') = Ueixi) + hC{xi) + 4'D(a;i) + 4'E(a:i) 

+h^O'l{x''), (3.7) 
iV0\x''))-' = Rf (xi) - hC.ix,) ~ a;^Di(xi) - x^F^ixi) 

+h^O^{x''), (3.8) 
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where 
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and 5^ : Q'' - 
functions which satisfy 



p 
, o, : (0,i) 





^ 0\ 




-pO 




i, 0/ 



1,2,3, O't : Q'^ 



b3x3 



(3.9) 
(3.10) 

(3.11) 
(3.12) 

(3.13) 

(3.14) 

(3.15) 

, fc = 1 , 2 are 



sup max |5''(x'')| < Co, 

0<h<ho a:''Gf2'' 



sup max ||o,f(xi)|| < Co, sup max ||(oJ')'(a;i)|| < Co 

0<h<ho^i<^(0,L) 0</i</io 

sup max max ||0^ ^ (x'')l| < Co, fc = l,2, 

0<h<ho 'J xi^efii^ 
for some constant Cq > 0. 
Proof. It can be easily seen 

t'\xi) = ei + he'2{xi)e2 + he':,ixi)e3-^{{e'2f + {9':,f)e, + h^o1ix"), (3.16) 

where ||o4||c2(o_l) < C. The relations (13. 5p and p.6p are the direct conse- 
quences of the relation p. 161) . Let us by u'^ : (0,i) denote the function 

u'^ = {l,h9'2,h9',f. (3.17) 

It is easy to see 

ye^x'^) = iu\x,) + x^^in'^Yixi) + xUb'^Yix,) \ n\x,) I b'^ix,)). (3.18) 
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The relations p.3p . p.7p . (j3.9L (I3.10p are the direct consequences of the rela- 
tions (j3.41 - (l3.6p and (j3.18D . The relation (j3.8p is the direct consequence of the 
fact that, for a regular matrix A and arbitrary B, which satisfies ||A^^B|| < 1 
(II • II is the operational norm), the matrix A + B is invcrtible and 

||(A + B)-i-(A-i-A-iBA-i)|| < M^"'B||2||A||-i 



-IIA-iBI 
To end the proof observe that 

Ci = RJCRJ , Di = R^DR^, El = R^TerJ^. 



Remark 1 By a careful computation it can be seen that , O3 and O4 (defined 
in the relation p.l6p ) are dominantly in 62,63 plane i.e. that we have for 
i = 2,3,4: 

o1{x,) = /f (a;i)62 + /f (xi)e3 + /ir'/(xi), (3.19) 
where ff, ff e 0^(0, L), supo<^<;,^^ ||r,''||ci(o,L) < C, for some C > 0. 
Remark 2 By a further inspection it can be seen that 



fl 




+ i03f), /I- 




(3.20) 




= P2{f2{0',9"] 




^3)')) -^2(^2^2 +^3P3), 


(3.21) 


/I 


= P3{f2{0',9"] 


)~lii92f + {^ 


^3)')) -^3(^2^2 +^^P3), 


(3.22) 


/I 




m-P3h{e'. 


0"), 


(3.23) 


/I 


= -P2ii9'2f + 


i9'3f)+P2f2i6 




(3.24) 



where f2{0',0") £ C^{0,L) is the expression that includes 0',9": 

f2{e\o") = \({P2+P3m'2f + {o'3?) + m + 9'3W2P2 + e',p,) 

-^{0'i)' + m?{0'2P2 + e',p:,f) 

Remark 3 It is not necessary to impose the condition (|2.9p . All we need is 
the existence of the expansions given by p.4p - p.6p . where P2,P3 G C^(0, L), 
including the statement of Remark [1] 

Remark 4 Although &^ makes the small perturbation of the central line, 
(xi, 0, 0), for xi G [0, L], it is not true that V©'' is close to the identity (like in 
the shallow shell model, see [31]). In fact, there are torsional effects of order 
on every cross section. This is the main reason why is the change of coordinates 
introduced in the next chapter useful. 
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4 r'-convergence 

We shall need the following theorem which can be found in |12) . 

Theorem 2 (on geometric rigidity) Let U C M™ be a bounded Lipschitz 
domain, m > 2. Then there exists a constant C(U) with the following property: 
for every v G VF^'^(J7; M™) there is an associated rotation R € SO(m) such 
that 

\\Vv - R||l2(c/) < CiU)\\ dist(Vt;, S0(m)||i2(c,). (4.1) 

The constant C(U) can be chosen uniformly for a family of domains which are 
Bilipschitz equivalent with controlled Lipschitz constants. The constant C{U) 
is invariant under dilatations. 

The following version of the Korn's inequality is needed. 

Lemma 1 Let C with Lipschitz boundary and u e L^iuj-^R"^). Let us by 
eij{u) denote eij{u) — ^{diU + dju). Let us suppose that for every i,j = 1,2 
we have that eij{u) G L'^{uj). Then we have that u G VF"'^'^(a;;R^). Also there 
exists constant C{uj), depending only on the domain uj, such that we have 

\\u\\w^.^{lj-M.^) < Ci^)i\ / udxidX2\ + \ / {xiU2 - X2Ui)dxidX2\ 

<t UJ J UJ 

+ E lleyNlli^M)- (4.2) 

Let us suppose that the domains lOs are changing in the sense that they are 
equal to uJs — A.sLU, where Ag £ M^^^, and there exists a constant C such that 
||As||, IIA^"'^!! < C. Then the constant in the inequality J^.iip can be chosen 
independently of s. 

Proof. The first part of the lemma (the fixed domain) is a version of the 
Korn's inequality (see e.g. ^^). The last part we shall prove by a contradiction. 
Let us suppose the contrary that for each n £ N there exists s" and u." G 
W^''^{uJs"',R?) such that we have 

I / tt"(ia;i(ia;2| + I / {xiu^ — X2U^)dxidx2\ 

+ < (4.3) 

where we have by ef^ (•) denoted the symmetrized gradient on the domain oj^n . 
Without any loss of generality we can suppose that 2(^^,1 ;R2) = 1- Let 

us take the subsequence of (s") (still denoted by (s")) such that Agn A 
and A-} A-i in ]R2x2^ 

Let us look the sequence m" = it" o A^n o A^^. It is clear that there exist 
Ci , C2 > such that 

Ci < ||<|1k.i.2(„^;K2) < C2, (4.4) 
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where we have put Woo ■— Aw. Thus there exists u G W^''^{uJoo',^'^) such 
that w" ^ u weakly in W^''^{uJoo',^'^)- Specially, by the compactness of the 
embedding ^ W^'"^ (see e.g. [2]), we also conclude the strong convergence 
It" — > M in L^(a;oo; R^)- Since it is valid Ag^ — ^ I, it can be easily seen that, 
from the weak convergence, it follows efj (it") o (A^™ o A~^) ^ e'^{u), weakly 
in i^(a;oo; K^), where we have by e^(-) denoted the symmetrized gradient on 
the domain Wqo- From the weak convergence we can conclude that 

||e,7(M)|U.(.^) <liminf||ej;(M")o(A,.oA-i)|U.(„^)=0, (4.5) 
for every i,j = 1,2. We can also from (14.31) conclude that 

/ udxidx2 =0, / {xiU2 — X2U\)dx\dx2 = 0. (4-6) 
Applying the standard Korn's inequality on the domain Woo, i.e. 



E 

tj=1.2 



||ey(M) - ey«)||L2(^^.R2)^, 



we conclude that m" -> u strongly in W^''^{ll>oo;^'^)- But then (|4.4I) . (j4.5p . 
(j4.6p make a contradiction with the version of the Korn's inequality (14. 2p on 
the domain Woo- 

Remark 5 The same proof can be done under the assumption that Us = Fs{uj), 
where Fs is the family of Bilipschitz mappings whose Bilipschitz constants we 
can control (i.e. the Lipschitz constants of Fg and are bounded by a univer- 
sal constant), provided that the family Fg is strongly compact in VF^'°°(w; R^). 
It would require more analysis to conclude the same result only for Bilipschitz 
mappings whose Bilipschitz constants we can control. 

Let us by a;' : M'^ -> denote the change of coordinates 

{x[,x'2,x'^) ^ x'{xi,X2,X3) ■.= 'Re{xi) \ ^2 \ ■ (4.7) 

By i?' we denote x'{f]) and uj'{xi) C R^ denotes x'{{xi} x w), for xi G [0,L]. 
The generic point in f2' is denoted with x' = (x'^, a;2, X3). Let us observe that 
by ([22]) and ((23t 



x'2dx'2dx';j — / X3(ix2(ix3 = 0, (4-8) 

'{xi) Ju'{xi) 

/ XjXgdxjdxg = / {x2 — x1)dx2dx3, (4-9) 

J u' (xi) J u 

h{lu) = {xj + xl)dx2dx3 ^ / {{x^)^ + {x'^)^)dx2dx'^, (4.10) 
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for all xi e [0, L\. By {d^y^) oO'^o P'') we denote diyj evaluated at the point 

In the sequel we suppose ho > 1 (see Theorem [1]) . If this was not the 
case, what follows could be easily adapted. Using theorem [5] we can prove the 
following theorem 

Theorem 3 Let G W^''^{t2'^;M?) and let 

E'' = ^[ dist^{Vy'\S0(3))dx. 

Let us suppose that 

limsup — ^ < +0O. (4-11) 

Then there exist maps R'' : [0,L] ^ S0(3) and R'' : [0,i] M^xs^ ^^^^^ 
|R| < C, R e tyi'2([o,i],M3x3) constants r'' G S0(3), c'' G smc/i 
</ia< f/ie functions y ' := (R )"^y'' — satisfy 

||(Vy'^) o o P'' - R"!U.(,,) < Ch\ (4.12) 

IIR'' - R''|U2([o,L]) < C/i^ ||(R'')'||i.([o,L]) < Ch, (4.13) 
I|R''-I||l~([o.l1) <C/i. (4.14) 

Moreover if we define 

^j, ^ ^ 1^ SNe^iL^l^,,,,,,, t = 2,3, (4.16) 
l_ 1 /'li(S3°e"''^''')-l3(52°6"oP'') 



i/ien, io subsequences, the following properties are satisfied 

(a) u'' in W^'^iO^L); 

(b) Vk in W^i'2(0, L), where Vk G W^'^{0, L) for fc = 2, 3. 

(c) ^ w weakly in W^''^{Q,L); 

(d) (Vv'^)°0"°P"-i ^ L^{Q), where A G W^'^{{),L) is given by 

A = i v'2 \ . (4.18) 

\v'3 w 

("ej sym ^^7^ uniformly on (0,L); 
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(f) the sequence "/^ defined by 



+4((«2 )' + e'2){^l) + + ^3)(^l) j , 

1 / {y'loO'^ oph){x)-hek~hx'^ 



h 

-vl{xi) - ) , for A: = 2,3, 



where {x')^ :— (0,— xjj,^^, is weakly convergent in L'^{Q) to a function 
7 belonging to the space C, where 

C = {7 e L^f2;R^) : / 7 = 0, 927,^37 e L^n;R^), 

J CO 

{^3l2{^i^ •) - X2'j3{xi,-))dx2dx3 = 0, for a.e. xi £ (0, L)}.(4.19) 



Moreover dkj'^ dkj in L'^{Q) for k — 2,3, 

Proof. We follow the proof of Theorem 2.2 in [25]. Applying Theoreni[2]as in 
the compactness result of [24] (using the boundedness of V6>'' and (V©''')^^ 
we can find a sequence of piecewise constant maps R'' : [0,L] — > S0(3) such 
that 

f ||(Vy'') o 0'' o P'' - R'^W^dx < Ch^, (4.20) 



and 



||R''(a;i + e) - R(a;i)f dxi < Ch'm + /^)^ (4.21) 



where /' is any open interval in (0, L) and ^ £ R satisfies |^| < dist(/', {0, L}). 
Let T] e C||f (0, 1) be such that 77 > and r^{s)ds = 1. We set r^h = Jiviji) 
and we define 

R''(a;i) := / 77ft(s)R''(a;i - s)ds, 

where we have extended R'' outside [0, L] by taking R''(xi) = R''(0) for every 
xi < 0, R'*(a;i) = R''(i) for every xi > L. Clearly ||R''|| < C for every h while 
properties ()4.13p follow from properties (I4.2ip . Moreover since by construction 
(see [H]) 

||R''(a;i + s) - R''(a;i)f < / dist^(Vy^ S0(3)) < C/^^ 
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for every \s\ < h we have by Jensen's inequality that 

||R'^-R"||i^([o^^]^„3x3 <C/i3 (4.22) 

By the Sobolev-Poincare inequaUty and the second inequaUty in (|4.13l) , there 
exist constants Q'* e R^^^ such that ||R'' - Q''||loo([o,2,]:R3x3) < Ch. Combin- 
ing this inequaUty with (|4.22p . we have that jjR'' — Q'*||l°°([o,l];K3x3) < Ch. 
This impUes that dist(Q'*, S0(3)) < Ch; thus, we may assume that Q'' belongs 
to S0(3) and by modifying Q'* by order h, if needed. Now choosing R — Q'' 
and replacing R'* by (Q^)^R^ and R.'' by (Q'')^R'', we obtain gJH). By 
suitable choice of constants c'* 6 M'^ we may assume that 

[ (y5'o6)''oP''-a;i) = 0, [ (y^ o©'' oP'* - /le^fc) 0, for fc = 2, 3. (4.23) 

Let A'' = By (HJil) there exists A G L°°((0, i); R^^^^) such that, up to 

subsequences, 

A'' ^ A weakly * in L°°{{Q,L);W^'"'^). (4.24) 
On the other hand it follows from (I4.13P and (|4.14D that 

A weakly in ^^^'^((O, L); R^^'''). (4.25) 

In particular, A e Vl^^'^((0, i); R^^^) and /i^^(R'' - I) also converges uni- 
formly. Using (j4.22p we deduce that 

a'' A uniformly. (4.26) 

In view of (j4.12l) this clearly implies the convergence property in (d). Since 
R'' e S0(3) we have 

A'' + (A'')^ = -hA'\A''f. 

Hence, A + A"^ = 0. Moreover, after division by 2h we obtain property (e) by 
(|4.26p . For adapting the proof to the proof of Theorem 2.2 in [3S] it is essential 
to see 

(Vy'') o &^ o P'' = (V(y'' o o F'')((V0'')-i o p'') 

= V/,(y ' o o P'')((V0'')-i o P''). (4.27) 

From KT7]\ it follows 

((Vy'') o 0'' o P'')((V6>'') o P^) ^ Vft(y'' o o p''). (4.28) 

and 

(yy^) oe^oP'^ -1= {Vhiy'" o&^oP'^ -0^o P'')){{V&^)-^ o P'' - Rf ) 
+{Vh{y' o o P'' - o P''))Rf . (4.29) 
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Let US notice that from (|2.8p . (|3.5p . (|3.6I) we can conclude 

6>fe = hOk + hx'k + Ofc(/i^) for fc = 2, 3, (4.30) 

where |10fe(/i3)||ci(r2) < C/i^. 

By multiplying (d) with (V0'')oP'» = Vft(6>''oP'') and using ([311), Km 
we obtain 

— ARe in ^2(12). (4.31) 

Property (b) immediately from (|43T|) by using ((O) . (|423)) and (|430l) . More- 
over, = Afei for fc = 2,3 so that Vk G W^2,2(q^^) g^j^^^g ^ e W^i'2(0,L). By 
using (e), (jXTOl) . (pi^ and (|OT|) from we conclude that 



^ sym(V„(y'' o o - o P^)RJ: 



< C (4.32) 

L2((0X);R3x3) 

The weak convergence of follows from p.7p . (|4.32p and the definition of Re. 
By using the convergence (j4.3ip and Poincare inequality on each cut {xi} x uj 
we can conclude 

/^^yj^^°®^°^ ^ 

^ (ARe)22.T2 + (ARe)23a;3 m L^H) . (4.33) 
By using (|23)) . (|4Jl) and (|4?30l) we conclude from (|433)) 
, _ {y'^ o ©'^ o P" - x'^ 1 



2 



h K^A 



f o e'^ o P'^ ^ A23,x'^ in L^{Q). (4.34) 



Let us note that since the left hand side of (|4.33l) i.e. (|4.34p is bounded 
in M^^'^(i7) the convergence in (I4.34p is in fact weak in W^^'^(f2). The only 
nontrivial thing to prove is the boundedness of diW2 in L^{f2). By the chain 
rule we have for i = 1,2,3 

d^{y\ o o p") = {{d^y\) o o p")((ai0i') o p") 

+((92^-) o ©'^ o P'')((ai02^) o P'') + ((aayf) o 0^ o P'')((5i03^) o p'')(4.35) 
and for fc = 2, 3 



d^{gl o 0" o p'^) = ;.[((9iy,f) o 0'' o p'')((afe0j) o p^) 

+ {{d2y-) o o P")((afc0^) o P") + ((agyj') o 0'' o P'')((afe0j) o P'')] (4.36) 

From and (|i35|) we conclude that the boundedness of diw!^ 

in L'^{Q) is equivalent to the boundedness of 

R^\ai0i + R52ai02 + R535i03 - /^(Paa^a - p'^x^) 



^2^(R2i/ ai0i+R2^y d,02 + BJi^j 9x03), (4.37) 
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in L'^{Q). By using (lO)) and (IXTSl) we conclude 

" h 



(4.38) 



The boundedness of Z2 in L'^{f2) is the consequence of (I3.5p . (I3.6P and (j4.14p . 
Now we have proved W2 ^ A23a:3 weakly in W^•^{^2). 
Analogously we conclude 



f y^^oe'^oP'^ ^-A23x'2,. (4.39) 



weakly in W'^''^{f2). Now, since w'' can be written as 

w\xi) = I (x>J - x'^w!^)dx2dx3, (4.40) 

it is clear that w'' converges weakly to the function w = — A23 = A32 in 
W^i'2(0, L). Let us define for /3'' : f]' ->■ , (3'' = j o x'-^ . By the chain rule 
we have 

d,f3l = (an'/) o {x')-' + {pW2+PW3)id27'n o {x')-' 

+ i-pW2+PW3)id3l'l)o{xr\ 

52/3," =^2(927-) ° (a^')"' -P3{dsj1) o (a:')"\ 

53/3," =P3(a27') ° i^T' +P2{dsj'l) o (XT'- (4.41) 

By differentiating (3i with respect to xj., with k=2,3, we have 

d2f3l = ^92(y^ o 0'' o o {x')-') + liiv^Y + e'2), (4.42) 
d^fil = ^93(yi' o 0'^ o P'^ o {x')-') + \{{v'i)' + e',)). (4.43) 

Let us analyze only 92/3i- We have by (j3.18L (|4.36p and the chain rule 

. ((9iy^) o 0" o P^ o (a;0-i)(p2^' - P3b\) 
d2f3, -2 

, {{d2yi) o 0" o pf^ o {x')-'){P2n'i -P3b^2) 



{{dsyl) o 0" o P^ o (a^0-i)(p2nt -P3b^) 
+ ^((«2)' + f?2)- (4.44) 
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By using ([S^]), (gSl), (EIHl), dUlI), dUll), (^351) and the definition of we 



can conclude tliat for proving the boundedness of 92 /^i it is enough to prove 
the boundedness of in L'^{fi) where 

■= + = ^^^2 + . (4.45) 

The boundedness of in L°°(J7) is then the consequence of the property (e). 
In the same way we can prove the boundedness of dzl3\. Using the Poincare 
inequahty and the fact that /^/^^^^^ I3i\dx2, dx'^ = 0, we deduce that there exists 
a constant C > such that 

{I3\{x)fdx2dxz <C f [{d2l3\{x)f + {dal3\{x)f]dx2dxs 

uj' {x\) Juj'(x-i) 

for a.e. xi € (0, L) and for every h. Ahhough the constant C depends on the 
domain, since all domains are translations and rotations of the domain uj, the 
constant C can be chosen uniformly. Integrating both sides with respect to xi, 
we obtain that the sequence (^^) is bounded in LF'^fi') so, up to subsequences 
I3\ (3-^ and dkl3i dkf3i weakly in L'^{f2'), for = 2, 3. From the relations 
(|4.4ip it can be concluded that j'l and dhJi dk-fi weakly in L'^{f2), 

for k — 2,3, where 7 = /3 o x' . For the sequences (/32), if^s)^ we have by 
differentiation that for j, k — 2,3 

djUl - (^^5, (yL' o e'^ o P" o {x')-') - h5,k - hw\l - <5,,)(-l)'= j . (4.46) 

By using the chain rule we see that for k = 2,3, 

d2iy1 ° 0'^ ° P" ° i^T') - h(^iid,yl) o 0'' o P'^ o {x')-'){P2n'l - p.b'l) 

+ {{d2y1) o e'^ o P>^ o {x')-'){P2n'i - pA) 
+ ((aayfc) ° 0' ° P'' o {xT'){P2n'i ~ P3b^s) 
o ©'^ o P" o {x')-') = h({{d,yl) o 0'' o P'' o (a^')"')(P3n? +P2b'i 
+ {{d2y'l) o 0'^ o P" o (a;')"')(P3n^' 
+ ((aayfc) ° 0" ° o {x')-'){P7^nl+p2bl 



h\ 



Now we want to check that for j,k — 2,3 

eMfi'')--^l{d,(3i + dkf3';). (4.47) 

is bounded in In the similar way as for f3i (relations (I4.44p and (14. 45^ *) 

we can using (13. 5p . p.6p . (j4.12D . (I4.14p and the property (e) conclude that for 
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every j, fc = 2,3, ejk{^^) G By using Korn's inequality (Lenima[T]) we 

have that there exists C > such that 

11/32 11^1-2(^^/(2:1)) + ll/^a l!vKi.2(tj'(2;i)) ^ C'f I / [32dx'2dx'r^\ 

+ 1 / f3ldx'2dx'^\ + I [ {x',(3'^ - x'2(3'^)dx'2dx',\ 

+ E l|e.fe(/3'')llL2(.'(.,)))> (4-48) 

for a.e. xi G (0, L). From the definition of and we see that the functions 
{I32{xi, •), I3^{xi, •)) belong to the space 

B^, = {/3 = (/32,/33) e W^'''{J{xi);^^) : / /Sd^dx^, = 0, 

[ (4/33 - x',f32)dx'2dx', = 0} (4.49) 

J uj' {x\) 

for every xi. By integrating (|4.48p with respect to xi we conclude that 132,(^3 
are bounded in i^(]7') as well as their derivatives with respect to X2,X3. From 
this we can conclude the same fact about 72 1 73 • The fact that the weak limit 
belongs to the space C can be concluded from the fact that for every h and 
a.e. xi {I32{xi,-),f33{xi,-)) e B^,. This finishes the proof of (f). 



4.1 Lower bound 

Lemma 2 Let , , E'\ R'', u'', , , 'y'\ (3'' = 7'' o {x')-\ 7, /3 = 
70 [x')^^, A be as in Theorem\^and let us suppose that the condition i4.11^ 
is satisfied and that 7'' ^ 7, 82'^'^ ^27, 937'' ^ 937 weakly in L'^{S7) i.e. 
f3'' 13, 52/3'' ^ 92/3, dsP'' 83/3 weakly in L^{n'). Let us define 

.?(44( ^^^-^°^^°^;i^-^-^^°^^^ -^-(^.) 

+4(«2')'(^i)+4(«3)'(^i)), (4-50) 
1 ({tkO&''oP'-){x)-&loph 



rilix) = 



/l2 



-D^(xi)-/i(4)^^"(a;i) , fork = 2,3, (4.51) 



and k'^ = T]'^ o (a;') ^. Then we have that r)^ r] weakly in L'^{fl) and 
dkV^ ^ dkV weakly in L'^{fl) i.e. k, 82^'^ 82^,, 831^''^ 83K weakly 
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m {[}'). Here 



»7i = 7i, 



V3 = 73 + f3^2 + fsX-i = 73 + .93-^2 + 532^3: 

K — r) o (x')^^ , 



(4.52) 
(4.53) 
(4.54) 
(4.55) 



fl are defined in Remark\^ and gj, can be easily defined for the above identities 
to he valid i.e. for k = 2,3, we define 



9k = P2fk - Psfk > 9k = Psfk + P2fk ■ 
The following strain convergence is valid 

and the symmetric part of G denoted by G, satisfies 

1 



G = syni(J-^A2 + K), 



where 



'J2"2 "T '^3^3 

wO'3 



1721/2 c^2''3 



K 



^2^2 -"^s^s 



X2W' 



82 K 



83 K 



(4.56) 
(4.57) 

(4.58) 

(4.59) 
(4.60) 



Moreover, 



liminf-ir / W^(x,Vy'')dx = \imM^ [ W{x^\{Vy^) o O''- o P^)dx'' 

^3{x, G{x))dx, 



1 



where Q3 is twice the quadratic form of linearized elasticity, i.e., 



Q3{x,F) 



9F2 



.(I)[F,F] 



(4.61) 



Proof. We follow the proof of Lemma 2.3 in [25 . Firstly, using Remark [1] it 
can be seen that 

Vi=li+hoi, 

V2=l2+ flx2 + /2a^3 + h02, 
V'3 = 73 + /3'2:2 + fix3 + h03, 
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where ||oi||(7i(i-2) < C, for some C > 0. The convergence of rj'* is an easy 
consequence of the convergence of 7''. The estimate (|4.12p imphes that the 
norm of G'* is bounded; therefore up to subsequences, there exists G G 
L^(J?; R'^^'^) such that (|4.57p is satisfied. In order to identify the symmetric 
part of G we decompose R'^G'* as foUows: 



(Vy')o0''op'^-i R^-i 



so that 



F'^ :^ sym ^ ° ^ ° ' ^ ^ sym(R''G'') + sym (4.63) 

The right hand side converges weakly to G + ^ by (|4.14p , (|4.57p and property 
(e) of the Theorem[31 Therefore the sequence F'' has a weak Umit in L^{0, L), 
satisfying F = G + To conclude we only need to identify F. Consider the 
functions 

From property (f) of Theorem[2]it follows that the functions M''+X2((f 2 )'+ 
^2) + ^sii'^^aY + ^3)1 which are equal to hj^ converge strongly to in L^{f2). 
Thus by property (a) and (b) of Theorem [3] we conclude that 

(jj'l ^ u ^ x'^iv'^ + e'^) - 4(w3 + in L^{Q). (4.65) 

By using the chain rule, the property (d) of Theorem [3l (|3.5L (13. 6L (I4.30p we 
can conclude that 

ai0^-Fu-a:2(e^P2+0^P3)'+a:3(e2P3-e3P2)'-^'^(ai4+(?i)-z;;,(5i4 + 0^,), 

(4.66) 

weakly in L'^{Q). From (j4.65p and (j4.66p we conclude that 

u' - ai4(«^ + e'^) - x'^iv'i + e'i) - d.x'M + e'^) - x'^ + e'i) 
= Fn - x2{e'^P2 + e'^p^y + X3{e'^p3 - e',p2y 

-v'2{dix'2+d'^)-v'^{dix'-, + 9'-,). (4.67) 

After some calculation we obtain 

Fn - ,/ + + v'^e'^ - x'^v'i - x'^vl (4.68) 

To identify F12 we have to do some straight forward computations. By using 
the chain rule, p.Sp . p.6p . Remark [U property (d) of Theorem Q we can 
conclude 

^91(^2' o 0" o P^) + -1 (p2d2{y\ o o P'^) - p.,d^{y\ o 0'' o P'- 
2F?2 + ^ - ^^'(^3 + dix'^) + , (4.69) 
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where lini;i_j.o ||05'||i2(^2:K3x3-) = 0. On the other hand it can be easily seen 
that 

92/3? =P2527!'-P3a37? 

+ ^((^2)' + ^2)- (4.70) 
From (|434)) . gUl), (|4J0)) we conclude 

2F?2 = -191(^2^ ° 0' ° - - ;^((^'2)' + O2) 

= diw!^ + d2f3'l + w0':i + wdix'.i-O'l. (4.71) 
By using (|4.34|) we conclude that the right hand side of (|4.71|) converges in 

w~^-^{n) to 

di{-wx'^) + 82(^1 + w9'^ + wdix'^ = -x'^w' + w9'^ + 82^1, (4.72) 

since f3i = Ki. On the other hand we know that the left hand side of (|4.7ip 
converges strongly in L^(J7) to 2Fi2 and thus we can conclude 

Fr2 = ^{~xW + we'3 + d2Ki). (4.73) 

In the same way one can prove 

Fi3 = ^ixW - wO'^ + ^3^^l). (4.74) 

To identify F22 let us observe that by the chain rule, p.Sp . p.6p and the 
property (d) of Theorem ([3]) we have 

^ (P292(y2 ° 0'' ° - 02 ° P") - P353(y2 ° ° P' " 02 O P" 

F^2-«2^2 + 02^ (4.75) 

where lim^_i.o IjOj ||l2(^2:R3x3) = 0. On the other hand we can conclude 

= ^ {P2d2{y2 O O P" - 02 O P") - P3a3(y^ O O P" - 02 O P''[ 

(4.76) 

In the same way as before we conclude that 

F22 -l'2^2+52K2. (4.77) 
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Analogously we can conclude 

F33 ^fs^s + Ss/^a- (4.78) 

To identify F23 = F32 we, by using the chain rule, p.5p . (|3.6I) and the property 
(d) of Theorem can conclude: 

^ (P292(y3 ° ° - 02 O P") - P3a2(y5 ° 0'' ° - 02 ° P")) = 

l(92y!,')o0''oP''~^;^0;, + O3\ (4.79) 
where lini^_i.o jjOg ||^2(^2:R3x3) = 0. In the same way we conclude 

^(P3a2(y2 ° 0" ° - 02 ° P") +P2a3(y5 O ©'^ O - 02 O P") 

^id2y^,)o0'^oP'^-v',9'^ + Ol (4.80) 
where lim/i^o ||C)4 ||l2(j-2;R3x3) = 0. It can be also concluded 

P2d2V3 - P3d3V3 = ^ {P2d2{y1 O &^ O P" - ©2 ° P'*) 

-V3d2{yl O 0" O P'' - 02 O P'^)) + (4.81) 
P352T72 + ^29377^ = ^ {pzd2{y2 ° ° P" - 02 ° P'*) 

+P2a2(yJ o 0" o P'' - 02 o P'')) - Jju;'\ (4.82) 

By summing the relations (|4.79p - (l4.82l) and letting ft, — J> it can be concluded 
that 

2F23 = t;2^^3+«3^2 + c'2K3 + c'3/«2. (4.83) 

To prove the lower bound we can continue in the same way as in the proof 
of Lemma 2.3 in [25^, by using the Taylor expansion, the cutting and Scorza- 
Dragoni theorem. 



4.2 Upper bound 

Theorem 4 (optimality of lower bound) Let u,w G W^'^{0, L) and Vk e 
W^'^{0, L) for k = 2,3. Let ■j be a function in C where 

C = {7 e L^n;R^) : f j = 0, 027,83^ € L^n;R^), 

(x372(xi, •) - a;273(a;i, ■))dx2dx3 = 0, Va;i £ (0, L)}. (4.84) 



Set 

G = sym(J - ^A^ + K). (4.85) 
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Here A, J, K are defined by the expressions ^.18^ , ^.59\ l and ^.60^ and rj, 
K are defined by the expressions ^.5^ - [4-55\ ). 

Then there exists a sequence {if'') C W^''^{(2'^ ,M.^) such that for u^,v^,w 
defined by the expressions ^.15\ l- (4-^^ the properties (a)-(d) of Theorem\3l 
are valid. Also we have that the property (f) of Theorem\^ is valid (which 
is equivalent that for rf^ defined by the expressions |^.5Q[ )-( f^.5j[ ) it is valid 
rj^ ^ r] weakly in L'^{Q) and dkTi^ dkr] weakly in L'^{fi) ). Also the following 
convergence is valid 



, J_ f w\x\Vy'')dx'^ Q^{x,G{x))dz 



lim — 

h 



tx (4.86) 

'fi'' ^ J n 

Proof. Let us first assume that UjWjVkjil are smooth. Then we define for 

(h'^u{xi) 
hV2{x,) 
hv3{xi) 

'X2{V2P2 + V'^p3){xi) - X'i{v'^P2 - v'^p3){xiy 

-h^ \ -X2{p3w){xi) - xz{p2w){xi) 

X2{P2w){xi) - Xz{pzw){xi) 



+h^rjix,,f,f), 



(4.87) 



where rj : f2 is going to be chosen later. The convergence (a)-(d) and 

that r]^ r] weakly in L?{Q) and dkrf^ dhf] weakly in L^(f?) can easily 
seen to be valid for this sequence. We also have 

' h'^u' ~h{v2P2 + v'^iPs) -h{v'3P2 - V2P3)' 



hv'2 



—hp^w —hp2W 
hp2W —hp^w 

I -X2{V2P2 + v'^PsY + X3{V2P3 - V3P2)' 
+ -X2{P3W)' - X3{P2W)' 

\ X2{p2Wy - Xsipsw)' 

+Oih^). 

From (|4.88p . by using (I3.8p . we conclude 





v'e' v'o' 



9211 



83 V 



(4.88) 





2 "3"2 ^3"3 , 

-X2{V2P2 + V3P3) + X3{V2P3 " W3P2) 
-X2{P3W') - X3{P2W') 
X2{P2W') - X3{P3W') 



82 K, 



83 K. 



(4.89) 
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Using the identity (I + M)^(I + M) = I + 2 sym M + M'^M we obtain 

{Vy'Yi^y'') = I + 2h^symJ + 2h^aymK + h'^A^ A + 0{h^), 

where HOC/i"^) ||l°°(j7;R3x3) < Ch^, for some C > 0. 
Taking the square root we obtain 

[(Vy'*)^(Vy'')]i/2 = I + h^G + 0{h^). (4.90) 

We have det(Vy'') > for sufhciently smah h. Hence by frame-indifference 
W{x, {Vy^) oB^o pf') = Wix, [Vy V{'^y^)Y''^ oO'^o P''); thus by (ICTl 
and Taylor expansion we obtain: 

^^W{x, (Vy'*) o o P^) ^ 1^3(0;, G(x)) a.e 

and by the property ii) of W for h smah enough 

^^W{x, (Vy'^) o 0^ o P>^) < ic(|| Jf + ||Kf + II Af) + Ch. 

The equahty (j4.86p follows by the dominated convergence theorem. Namely, 
we have 

/ W'^ix, Vy'')dx =^ f W{x, (Vy'') o ©'^ o P^)dx 
^ o / Qz{x,G)dx. 

In the general case, it is enough to smoothly approximate u, w in the strong 
topology of W^'^, Vk in the strong topology of W^''^, and t], dkrj in the strong 
topology of and to use the continuity of the right hand side of (|4.86p with 
respect to these convergences. 

Remark 6 Notice that 

K = I A' I a;'2 I | d2i^ \ | = L + ( A' ( xM I ^2/3 | ^3/3 



X' 



3 , 



Here f3 = -f o [x') ^ and 

/O 

L = gl gl I . (4.91) 

From the fact that 7 e C we can conclude (3 ^ B, where 

B = {l3e L^{n';R^) : j /3 = 0, 92/3,93/3 G L'^{Q';R?), 

•J CJ 

/ (a;^/32(a;i,-) -X2/33(xi,-)dx2dx3 =0, for a.e. xi 6 (0, i)}.(4.92) 

J u)'{x{) 
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4.3 Identification of tlie I^-limit 

Let Q : (0, L) X R X so(3) [0, +oo) be defined as 
Qixi,t,F) = 

Q3{x,{F{x'2]+te 



min / 



do a 



d^a dx^dx'^, 



(4.93) 

where Q3 is tlie quadratic form defined in (I4.6ip . For u,w E W^''^{0,L) and 
V2,V3 € W^'^{Q,L) we introduce the functional 

I°{u,V2,V3,w) := i l\{x,,u' + v'2e'2+v',e', + ^{{v'2)' + {v',f),diA)dxu 

(4.94) 

where A e ^^^'^((0, L); so(3)) is defined by ffTJ^ . We shall state the result 
of P-convergence of the functionals -^l'^ to /°. Before stating the theorem we 
analyze some properties of the limit density Q. 

Remark 7 By using the remarks in the beginning of chapter 4 in |25j the 
following facts can be concluded: 

a) The functional Qaix, G) is coercive on symmetric matrices i.e. there exists 
a constant C > 0, independent of x, such that Q3{x,G) > C||symGp, 
for every G (this is the direct consequence of the assumption iv) on W). 
The minimum in (|4.93|) is attained. Since the functional Qsix, G) depends 
only on the symmetric part of G, it is invariant under transformation a n> 
a + ci + C2{x')^ and hence the minimum can be computed on the subspace 

V^, := \a eW^-^{uj'{xi),R^) : [ a = 0, 

I (2:302 — 2; 2 03)^2; 2^2:3 = >. 

Ju>'{xi) J 

Strict convexity of (33(2:, •) on symmetric matrices ensures that the mini- 
mizer is unique in V . 

b) Fix xi £ {0,L),t eR and F £ so(3). Let o™*" G F be the unique minimizer 
of the problem (gjlS]). We set 







( 2;^ I +tei, 6^fc = ^-_^(a;,l), 

and we call B'^'^ the matrix in R'^'*^ whose elements are given by {B^'^)ij = 
b^j' . Then a™*" satisfies the following Euler-Lagrange equation: 

[ ^ (B'''=9fea™'",a^^)dx^dx:5 = - /" Y.(^''^9,dh^)dx'2dx'3, 

(4.95) 
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for every Lp G W^''^{uj' {xi);M.^^^). From this equatfon it is clear that a™*" 
depends linearly on {t,F). Moreover Q is uniformly positive definite, i.e. 

Q(a;i,i,F) > C(i2 + ||Ff ), Vi e M, VF e so(3), (4.96) 

and the constant C does not depend on xi. 

c) By mimicking the proof of Remark 4.3 in [3S] it can be seen that there 
exists a constant C (independent of xi, t and F) such that 

||92a™™||L2(^/(^j).R3x3) + ||93a"""|lL2(„,(^^).R3x3) < C"||5||^2(„,(^J.R3X3), 

(4.97) 

for a.e. xi € (0,L). To adapt the proof we only need to have that the 
constant in the Korn's inequality 

(4.98) 

can be chosen independently of xi. This is proved in Lemma [T] 

d) When Qa does not depend on X2, x^ we can find a more explicit repre- 
sentation for Q. More precisely Q can be decomposed into the sum of two 
quadratic forms 

Q(xi,i,F) = Qi(xi,t) + Q2(a:i,F), 

where 

Qi(xi,t):= min Q3(2:i, (tei|a|6)), (4.99) 
Q2(xi,0,F) :=g(xi,0,F). (4.100) 



The relations (|4.8p are only needed for this. If we assume the isotropic and 
homogenous case i.e. 

Q3(F) = 2m|^-!^|' + A(traceF)2, 

then after some calculation (see Remark 3.5 in 24 ) it can be shown that 
;3A + 2/ 

^^fF„ / 

u'{xi) 

2X^da 



^ M3A + 2^^2 
A + ^ 

Q2(a;i,F) = ^ii^A±M(Fi2 / {x',rdx',dx'. 



A + /i 

+2F12F13 / Xa^gdxjdxg + Fi3 / {x'^)'^dx'2dx';^ 

J uj' {x\ ) J bj' {x\ ) 

+MTr23, 

where the constant r is so-called torsional rigidity, defined as 

r(w'(xi)) = t(w) = / (a:2 + ^3 ^ X2d3(p + x^d2^p)dx2dx^ 
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and Lp is the torsion function i.e. the solution of the Neumann problem 

V93 = in w 

dyip = — (2:3, —X2) ■ V on duj 

The following theorem can be proved in the same way as Theorem 4.5 in [53] 
(we need Theorem [31 Lemma [21 Theorem [H Remark [5] and Remark [7]). 

Theorem 5 As h 0, the junctionals -^l'^ are F-convergent to the func- 
tional I*^ given in 1^4.94^ , in the following sense: 

i) (compactness and liminf inequality) if limsup^_j.Q ft.^'*/'' < +00 then 
there exists constants R G S0(3), c'^ G such that (up to subsequences) 
R — > R and the functions defined by 



r ■■= (R )^ - C\ -jj h^ -dX2dX3 

vu-^j^ h ^^"^^"3 

satisfy 

(a) (Vy^) oO'^ oP^ ^ \ m L^{f2). 

(b) there exist u,w ^ W^'^{0, L) such that u and w weakly in 
VKi'2(0,i). 

(c) there exists Vk G VK^'^(0,i) such that v'^ — > Vk strongly in W^'^{0,L) 
for fc = 2,3. 

Moreover we have 

liminf ^I'^iy'^) > I^'iu, V2,V3,w). (4.101) 

a) (limsup inequality) for every v,w £ W^'^{Q,L), V2,V3 £ W'^'^{0,L) there 
exists {y^) such that (a)-(c) hold (with y replaced by y'^j and 

hm ^I\y") = l''{u,V2,V3,w) (4.102) 
Remark 8 Let /2, /s G i^(0, L). Wc introduce the functional 

J° = I°{u,V2,V3,w)- f ^ /fci-fc, (4.103) 

for every u G W^'^{0,L), V2,V3 € W^^^{0,L), and w G W^'^{0,L). The func- 
tional J° can be obtained as /^-limit of the energies -^l'^ by adding a term 
describing transversal body forces of order h^ (see [T3], see also [32] )■ For lon- 
gitudinal body forces see [IH] . The problem for longitudinal body forces arises 
because the longitudinal forces should be of order /i^, the same order as for 
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the model in [M] . One needs to impose certain stability condition to see which 
model describes the behavior of the body for the longitudinal forces of order 
h\ 

Remark 9 The term u' + W2^'2 + ^^3^3 + 5(("2)^ + W^Y) the strain measures 
the extension of the central line (which is of the second order) . Namely, if we 
approximate the deformation of the weakly curved rod by: 

(pi{xi,X2,X'i) = xi + h^u + h^x'^{v'^ + e'^) + h^x'^{v'^ + 9'^) (4.104) 
(pkixi,X2,X3) = h6k + hx'^ + hvk + h^{x'i.)-^w, forfc = 2,3, (4.105) 

we see, that it is valid 

\\di(p{xi,0, 0)f - ||ai0''(a;i,O,O)f = h^(2u' + 2v'^0'^ + 2v'^e'^ 

Remark 10 The existence of the solution for the functional under the 
Dirichlet boundary condition for Vk at both ends of the rod can be proved 
directly. It is also enough that we impose W2 , ^2 , ^'s , at the one end. The 
existence can also be proved for the free boundary condition under the hy- 
pothesis that ff.dxi — 0, J^^ xif /.dxi = for k — 2,5. It can be done in 
the same way as the proof of Lemma 5 in [32] . 
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